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System  Complexity  translates  to  EM  Complexity 


■  EM  Complexity 

—  Geometric  complexity  and  distributed  nature  of  the 
packaging  hierarchy  (“coupling  path”) 

—  Broad  frequency  bandwidth  of  the  interfering  signal 

—  Non-linearity  of  the  terminations 

■  Tackling  Complexity 

—  Hierarchical  approach  to  the  modeling  of  electromagnetic 
interactions 

•  From  lumped  models,  to  transmission-line  models,  to  quasi-3D 
full-wave  models,  to  3D  full-wave  models 

■  Abstracting  Complexity 

—  Systematic  order  reduction  of  numerical  models  of  the 
coupling  path 

—  Equivalent  circuit  representation  of  the  coupling  path  for  its 
seamless  incorporation  in  network-oriented  non-linear  circuit 
simulators 
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SUMMARY  OF  MODELING  OBJECTIVE 
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Specific  Subtasks 


■  Task  2.1  Development  of  a  coupling  path  modeling 
methodology 

■  Task  2.2  Development  of  a  (EMI)  source  modeling 
methodology 

■  Task  2.3  Non-linear  Transient  Simulation  of  the 
hybrid  lumped-distributed  non-linear  network 
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Subtask  1:  Modeling  of  the  Coupling  Path 
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Subtask  1 :  Modeling  of  the  Coupling  Path  (Cont) 


■  Our  modeling  approach  is  hierarchical 

—  EM-Physics  driven  “divide  and  conquer”  approach 

■  Use  suitable  EM  modeling  approach  for  individual 
blocks 

—  RLC  models  for  short  interconnect  at  chip  &  package  level 
—  Multi-conductor  Transmission  Lines  (MTL) 

•  Balanced  Interconnects  at  the  board  level 

—  Quasi-3D  EM  Modeling  for  PCB  power  delivery  network 
—  Full-wave  modeling 

•  Unbalanced  interconnects  at  the  board  level 

•  Coupling  through  cables 

•  RF/microwave  packages  &  boards 

■  Reduction  of  EM  models  to  non-linear  network 
simulator  (e.g.  SPICE) 

—  Direct  model  order  reduction 
—  Equivalent  circuit  synthesis 


■ 
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Block  Diagram  of  EM  Modeling  Flow 


Quasi-3D  Modeling  of  Power  Delivery  Network 


EM  field  between  power/ground  plane 
pairs  exhibits,  approximately,  a  two- 
dimensional  variation 
•  Use  simple  2D  FDTD 
Three-dimensional  features  (vias,  pins, 
slots,  voids,  etc)  require  locally  3D 
models  to  capture  the  correct  physics 


jwsE^  (i,  j)  =  Ax-'  {i  + 1,  j)  -  (i,  j)) 
-jw^AyH^  (i,  j)  =  £  (i,  j)  -  £  {i,  j  - 1) 


jco^AxH^  (i,  j)  =  £  (i,  j)  -  £  {i  - 1,  j) 
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The  resulting  discrete  EM  model 
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y(s)  =  F^{A  +  sBy^FV(s)  where  A  = 


B  = 


rri  /  _ \  ”  1  _ 

^he  multiport  transfer  function  is:  ¥(5)  =  F  A  +  j  F 


•  Direct  time-domain  simulation  is  possible  if  desired 

•  Could  be  synchronized  with  time-domain  IE  solver 

•  Alternatively,  model  order  reduction  of  the  transfer  function 
using  an  Arnoldi-based  subspace  iteration  method  (PRIMA) 
is  used  for  the  generation  of  a  compact  multiport 
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Multi-port  representation  of  the  reduced-order 
model  is  in  terms  of  rational  functions 


■  Form  compatible  with  circuit  simulators  that  support 
rational  function  models  (e.g.  H-Spice,  ADS) 

■  Equivalent  circuit  synthesis  is  possible  also 

■ 
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Passive  synthesis  through  Foster  forms 


For  a  passive,  reciprocal  multi-port,  its  admittance  (or  impedance) 
matrix  may  be  cast  in  the  following  form: 

a  a  ^ 

ll(s)  =  G,+sC,  +  ^  + 


=l[S-Pa  ^-Pc 


where  it  is: 


Re{jt7j<0,  q  =  \,2,....Q 

matrices  and  ^  =  0, 1, 2, are  real  and 

positive  semi-definite 

Re{aj>0,  q  =  \,2,....Q 

0<lm{a^ }  Im{/7 J  <  |Re  ^ }  Re  {/? j| ,  ^  =  1, 2, . . .  .g 
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Equivalent  circuit  for  a  two-port 
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Two-port  sub-circuit  for  use  in  H-SPICE 
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The  advantages  of  the  direct  use  of  the  rational 
function  macro-model  in  the  simulator 


H-SPICE 


Equivalent 


3  ports 

1 0  ports 

36  ports 

0.7  sec 

5.73  sec 

21 .2  min 

2.3  sec 

676.55  sec 

>12  hours 

H-SPICE 


Equivalent 


3  ports 

1 0  ports 

36  ports 

25.9  KB 

53.4  KB 

172  KB 

308.2  KB 

269  MB 

>2  GB 

The  secret  is  in  the  recursive  convolution  that  is  utilized  for  the 
interfacing  of  the  macro-model  with  the  rest  of  the  circuit 
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Application  to  a  4-layer  PCB 


8  cm  X  4  cm,  4  layer  board.  A  2mm  gap  is  present  all  the  way  across  the 
top  ground  plane.  A  total  of  58  pins  are  used. 


Top  view 


Cross-sectional  view 
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|Z21|  From  H-SPICE 
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Application  to  a  4-layer  PCB  (cont’d) 


H-SPICE  calculation  of  the  transient  response. 

Switching  occurs  at  port  1  while  port  2  is  terminated  with  50ohm. 
Voltage  at  port  1  is  depicted  by  the  solid  line. 

Voltage  at  port  2  is  depicted  by  the  dotted  line. 
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Validation  study  in  progress 


Intel  4-laver  test  board 

1®‘  layer:  Power  plane 
2"'*  layer:  Ground  plane 
S’’'*  layer:  Ground  plane 
4*'^  layer:  Power  plane 
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Samples  of  the  generated  mesh 
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Generated  multi-port  for  switching  noise 
simulation  in  SPICE 
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Rational  function  multi-port  &  equivalent  circuit 
synthesis  from  raw  data 


■  Data  either  from  measurements  or  EM  field 
solvers 

■  Step  1 :  Rational  function  fitting 

—  Process  guarantees  stability  but  not  passivity 
—  Check  fitted  multi-port  for  passivity 

•  If  not  passive,  constrain  fitting  using  Foster  constraints 

•  Repeat  fitting 
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Validation  of  rational  function  fitting 


PCB  interconnect  test 
structures  courtesy  of  Intel 
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Validation  of  fitting  (cont’d) 


Measured  |YgJ*  Synthesized  |YgJ 


(*)  Measurements  courtesy  of  Intel  Corporation 
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Task  2.1  Summary 


■  EM  modeling  flow  for  the  coupling  path  established 

■  Quasi-3D  EM  Model  for  power  delivery  network 
completed 

—  Validation  in  progress 
—  Further  enhancements  include: 

•  Incorporation  of  hooks  for  balanced  MTLs 

•  Incorporation  of  hooks  for  matrix  transfer  functions  in  Foster  form 

■  Capability  for  SPICE-compatible  broadband  multi-port 
macro-model  generation  in  place 
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Subtask  2.3:  Non-linear  Transient  Simulation 


■  Physics-oriented  nonlinear  transient  simulator 

■  Full  wave  modeling  of  printed  circuit  boards  (PCBs)  with  fine 
geometric  features,  finite  (and  possibly  inhomogeneous) 
dielectrics,  and  nonlinear  loads  and  circuits 

■  Interfaces  with  SPICE  solvers  and  models 


June  2002 


Introduction  (cont.) 


Time  Domain  Integral  Equation  (TDIE)  based 
Marching-on-in-time  (MOT)  solvers 


Have  been  known  to  the  acoustics  and 
electromagnetics  communities 
since  the  sixties. 


Compared  to  frequency  domain  solvers, 
they  can  solve  nonlinear  problems  directly 

Liu  and  Tesche  (1976), 

Landt  and  Miller  (1983), 

Djordjevic  and  Sarkar  (1985), 

Deiseroth  and  Singer  (1995), 

Orlandi  (1996) 
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Introduction  (cont.) 


However  TDIE  -  MOT  solvers  have  long  been  conceived  as 

■  Unstable:  ...  But  recently  many  proposals  have  surfaced  for 
stabilizing  these  schemes  (Davies,  Rao  and  Sarkar,  Walker, 
Smith,  Rynne,...); 

■  Siow:  . . .  Computational  complexity  has  prohibited 
application  to  analysis  of  iarge  scale  problems! 

. . .  PWTD  technology  removes  this 
computational  bottleneck 


0(N,  N,  log^  N,) 


Now,  we  can  /  should  be  able  to  rapidly  solve  large 
-scale,  nonlinear  problems  using  the  PWTD 
technology!!! 
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Formulation  -  Problem  Definition 


■  Given 


■  an  of  3-D  inhomogeneous  dielectric 
bodies  (F),  and  3-D  arbitrarily 
shaped  PEC  surfaces,  wires,  and 
surface-wire  junctions  (S ), 

■  Multiple  linear  /  nonlinear  lumped 
circuits  connected  to 

■  a  temporally  bandlimited  excitation 

■  Solve  for 

■  all  transient  currents  and  voltages 
induced  on  the  interconnect  geometry 
vSljF  and  the  the  circuits  (ckt  1,  ...,  ckt  m) 

■  radiated  electric  and/or  magnetic 
fields  if  required 


◄ -  - ► 


distributed  lumped 
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Formulation  -  Time  Domain  Electric  Field 


■  Assume  spatially-variable, 
frequency  independent 
permittivity  and  free-space 
permeabiiity  in  v,  and  thin 
wires  in  S . 

■  Radiated  eiectric  fieid: 


E''(r,r)  = - A{Y,t)  +  c^\dtyv  ■A{r,t') 


(1) 


0 


A(r,0  =  ^ 

Att 


I  ^  (r 0  +  \  dV  *  3^,,,  (r 0 


,  (2) 


R  = 


r-r 


'diel  (r.  0  =  (^(r)  -  ^0 )  T-  E'""'  (r,  t) 

ot 
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Formulation  -  MOT  Algorithm 


Ns  Nr 

©  Expand  current  as  J(r ,  0  ~  ^  ir)Tj  (t),  q  =  s,w,sw,d 

n=\  j=0  ,7-  ^  ^  ^ 

~  ^ PEC  ^DIEL 

■  Surface  and  Volume  Spatial  Basis  Functions 
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Formulation  -  MOT  Algorithm  (cont.) 


Construct  a  system  of  equations  by  applying  spatial  Galerkin  testing 
at  each  time  step. 


for  the  j'thWme  step: 


i-l 


Zol,-  =  V? 
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Formulation  -  Transient  Circuit  Simulator 


■  SPICE-like  transient  circuit  simulator 

■  Performs  linear  and  nonlinear  large-signal  analysis  using 
Modified  nodal  analysis 

■Nonlinear equations  solved  using  multi-dimensional 
Newton 

■  Incorporates  the  following  circuit  elements: 

■  Independent/dependent  voltage  and  current  sources 

■  Resistors,  inductors,  capacitors 

■  Diodes,  BJTs,  MOSFETs,  MESFETs 


— If 


o 
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Formulation  -  Transient  Circuit  Simulator  (cont.) 


Circuit  behavior  described  by  time-domain  nodal  equations 
(using  trapezoidal  rule  for  the  fih  time  step) 


immediate 

interactions 

unknowns 

sources 

^  ) 

previously 

computed 

interactions  from 
previous  time  step 

For  m  circuits,  the  nonlinear  system  of  equations  is 
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Formulation  -  Coupled  EM  /  Circuit  Equations 


■  Combine  the  EM  and  circuit  nodal  equations  into  a  single  consistent 
system  to  be  solved  for  EM  and  circuit  unknowns  at  each  time 

step 


■  Coupling  between  the  circute  and  th^e  interconnect  structure  is 
accounted  for  by  matrices  C  and  C 

V  i 
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Formulation  -  the  MLPWTD  Algorithm 


for  typical  PCB 
structures 


Level 

Direct(O) 

1 

Source  block 

□ 

Observer  block 

■ 

■ 
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Analysis  of  Active  Nonlinear  Microwave  Amplifier 


■  Objective 

Characterize  structure 
from  2-10  GHz 


Simulation  Parameters 

TV  =1468,  TV  =7096 

*3  K 

TV^+7V^=8564,  tv,  =512 
=  5.0  ps 

'  Excitation  pulse 

cDq  =  2;z'x6x10^  rad/s 
t'  =  t-6<j,  cr  =  6. 82xl0~^^s 
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Amplifier :  Nonlinear  circuitry 


■  MESFET  Circuit  Model 


4  (Vgs  >  V*  )  =  (4,  +  4  +  A^vl  +  4Vg^ )  tanh(av^^ ) 
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amplitude  (V) 


Amplifier:  Plane  Wave  Interference 


■  Transient  input/output  waveforms 


■  MOT  results* 


S21  -  MOT 
S11  -  MOT 
S21  -  REF 


O  S11  -  REF 


23456789  10 

Frequency  (GHz) 


*  C.  Kuo,  B.  Houshmand,  and  T.  Itoh,  “Full-wave  analysis  of  packaged  microwave  circuits  with  active  and  nonlinear 
devices:  an  FDTD  approach,”  IEEE  Trans.  Microwave  Theory  Tech.,  vol.  45,  pp.  819-826,  May  1997. 
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Amplifier  +  Shield:  geometry 


■  Objective 

Determine  effect  of  shielding  box 
on  S-parameters 

■  Geometry 


Simulation  Parameters 

=  3288,  =  7096 

A^^+A^^  =10384,  A^,  =700 
=  6.25  ps 


PEC  shieRnhg 
(640xl86x690mils) 
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Amplifier  +  Shield:  S-parameters 
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Amplifier  +  Shield: 


W“{Y,t)  =  E„  cos{co/)e~''"l^’^ 
E„=500x  +  500y  F/m 

<^0  =  2;7'  X  6  X 10^  rad/s 
t'  =  t  -6(j  +  (z-y)I  c 
cr  =  8.68xl0““s 
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Plane  Wave  Interference 


Transient  response  at  transistor 

output 
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Amplifier  +  Shield:  Plane  Wave  Interference 


E“‘(r,0  =  E„ 

=  2;t X 4 X lO’  rad/s 
^'  =  ^  — 6(j  — (k-r)/c  (J  =  l/4;Tns 

e,  =  145,  =  0 

k  =  0.5736x-0.8192y 


Transient  response  at  transistor 

output 
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Amplifier  +  Shield:  Plane  Wave  Interference 


E“"(r,0  =  Eo  cos{co/)e~'"l^'^ 

=  2;^  X  3  X 1 0^  rad/s 

('  =  (-6(7-(k.r)/c  cr  =  3/14;rns 

6;  =115,  ?i,  =45 
k  =  0.6409(x  +  y)  -  0.4226z 

Eg  =  -500  V/m,  Eg  =  500  V/m 


Transient  response  at  transistor 

output 
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Amplifier  +  Shield: 


E‘“(r,/)  =  E„  c,o&{co/)e‘''l^'’' 

o)(^  =271x6x10)^  rad/s 

t' —  t  —  6a —  ct  — 3/1  l;7r  ns 


6».  =180,  ^.=0 
k  =  -z 

=  -500  V/m,  =  500  V/m 
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Voltage  (V) 


Plane  Wave  Interference 


Transient  response  at  transistor 

output 


Task  2.3:  Summary 


1)  A  MOT  algorithm  based  on  a  hybrid  surface/volume  time  domain 
integral  equation  has  been  developed  for  analysis  of  conducting/ 
inhomogeneous  dielectric  bodies, 

2)  This  algorithm  has  been  accelerated  with  the  PWTD  technology 
that  rigorously  reduces  the  OiN^Nl)  computational  complexity  of 
the  MOT  solver  to  0(#y,#^log^Af^)^for  typical  PCB  structures, 

3)  Linear/Nonlinear  circuits  in  the  system  are  modeled  by  coupling 
modified  nodal  analysis  equations  of  circuits  to  MOT  equations, 

4)  A  nonlinear  Newton-based  solver  is  used  at  each  time  step  to 
consistently  solve  for  circuit  and  electromagnetic  unknowns. 

5)  The  proposed  method  can  find  extensive  use  in  EMC/EMI  and 
signal  integrity  analysis  of  PCB,  interconnect  and  packaging 
structures  with  realistic  complexity. 

■ 
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